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TRANSONIC FLOW PAST CONE CYLINDERS ‘ 

By George E. Solomon 


SUMMARY 

Experimental resulti< are prei^ented for tranmnic How ])a^st 
cone-cy Under, axially ayinmetric bodies. The drag coefficient 
and surface Mach number are studied as the jree'-stream Mach 
number is varied and, wherever possible, the experimental re- 
sults are compared with theoretical predictions. Interjero- 
metric results for several typical flow configurations are shown 
and an example of shock-free supersonic to subsonic compres- 
sion is experimentally demonstrated. 

The theoretical problem of transonic dow past finite cones 
is discussed briefly and an approximate solution of the axially 
symmetric transonic equations, valid for a semi-infinite cone, 
is presented. 

INTRODUCTION 

Transonic flow past certain two-dimensional bodies has 
been the subject of sev^eral recent papers and the phenomena 
are well understood. The theoretical results of Cole (ref. 1), 
Guderley and Yoshihara (ref. 2 ), Vincenti and Wagoner 
(ref. 3), and others apply to two-dimensional symmetrical 
double-wedge airfoils. The experimental results of Bryson 
(ref. 4) and Griffith (ref. 5) substantiate the theoretical work 
in a very satisfactory manner. More recently, Vincenti 
and Wagoner (ref. 6) and Guderley and Yoshihara (ref. 7) 
liave analyzed the transonic flow past two-dimensional un- 
synimetrical sections, that is, lifting double-wedge airfoils. 
Current experiments on lifting double-wedge airfoils (ref. 8) 
at the Guggenheim Aeronautical Laboratory of the California 
Institute of Technology indicate that agreement between 
theoretical and experimental results will again be obtained. 

Two-dimensional and axially symmetric transonic flows 
are of considerable theoretical and practical interest since 
these two specialized problems are limiting cases of the more 
complex problem of the flow about an arbitrary three- 
dimensional body. 

The study of axially symmetric transonic flow is not so 
complete as that of two-dimensional flow. In recent years 
several papers, notably those of Von Kannaii (ref. 9) and 
Oswatitscli and Berndt (ref. 10), have studied the similarity 
rules of axially symmetric transonic flovc. . Also, Yoshihara 
(ref. 11) has calculated the flow about a finite cone at a free- 
stream Mach number of 1 by a relaxation technique and has 
obtained some experimental verification of the theoretical 
result. The hodograph problem for general transonic; flow 
past finite cones is discussed in reference 8. However, theo- 
1 ‘etical solutions or experimental results for tlie complete tran- 
sonic regime ai*e not, at present, available. The present 


paper ])n'sents the results of an experimental investigation 
of the transonic flow past cone-cylinder bodies. A conical 
tip followed by a cylindrical afterbody was chosen as the 
experimental model for two primary reasons: (1) The 
relatively simple geometry of a cone-cylinder body may 
simplify the theoretical problem, and (2) viscous effects are 
minimized; that is, the boundary layer on tiie cone surface 
is in a region of decreasing or constant pressure so that the 
presence of the boundary layer will not greatly alter the shape 
of the body forward of the cone shoulder. 

Theoretical results 'for the supersonic flow past a cone 
were first presented in 1929 by Busemann (ref. 12). Buse- 
mann s solution postulates a semi-infinite cone and assumes 
that the flow is conical; that is, along rays thi’ough the apex 
of the semi-infinite cone, the flow parameters such as pres- 
sure and velocity are constant. The solution is found by a 
geometrical construction in the hodograph plane and it is 
readily apparent that a conical solution exists only so long 
as a shock wave is attached to the cone apex and, therefore, 
the free-stream Mach number is supersonic. It is interesting 
to note that Busemann’s solution predicts smooth shock-free 
compression from supersonic to subsonic flow for particular 
combinations of cone angle and free-stream Mach number. 
The conical solution also shows that for a given cone angle 
and free-stream Mach number Moo the surface Mach number 
is always less than the Mach number immediately behind 
the conical shock wave; as decreases, the surface Mach 
number decreases and eventually passes from supersonic to 
subsonic values. As was mentioned by Busemann, the 
conical solution for a semi-infinite cone is completely valid 
for a finite cone so long as the flow is everywhere supersonic, 
but when the surface Mach niunber is less than sonic the 
perturbation due to the corner or shoulder of the finite cone 
is propagated forward through the subsonic portion of the 
field destroying the conicity of the flow. Thus, the Buse- 
mann solution is completely valid for a finite cone only so 
long as Moo is large enough so that the surface Mach number 
is greater than sonic. 

Taylor and Maccoll (ref. 13) in 1933 presented the results 
of a numerical integration of the axially symmetric equations 
of motion for conical flow about semi-infinite cones and also 
presented experimental verification of their theoretical re- 
sults. Further experimental verification by Maccoll (ref. 14) 
was published in 1937. Both of the above papers noted that 
deviations of the experiments from the theoretical predic- 
tions, notably in the shape of the shock wave, are apparent 
when the surface Mach number is subsonic. 


» Supersedes NAG A TN 3213, “'J'ransonlc Flow Past Cone Cylinders” by George E. Solomon. 19M. 
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As was mentioned previously, Yoshihara (ref. 11) has 
computed, by relaxation methods, the flow about a cone 
cylinder at M^ = 1.00 and has experimentally verified the 
calculation. However, theoretical solutions do not exist for 
the complete transonic regime. Solutions have not been 
developed for the flow past a finite cone when is sub- 
sonic or when is between sonic and the value of at 
which Busemann’s conical solution becomes valid. Drougge 
(ref. 15) has computed the flow field between a detached 
shock wave and finite cone by relaxation methods; however, 
the position and shape of the detached shock wave were 
determined initially from schlieren photographs. 

The experimental results reported in this paper cover 
several interesting features of the transonic flow about finite 
cones. The deviations of the surface Mach number from 
the values predicted by conical theory are examined for 
values of such that the flow field is transonic in nature. 
The behavior of the surface Mach number for subsonic values 
of M^, and as approaches sonic from subsonic values, 
is examined in some detail so that an extrapolation to = 
1.00 may be made. The above surface Mach number data 
lead naturally to the evaluation of the drag coefficient, and 
experimental values of the drag coefficient in the transonic 
regime are presented. 

The phj-sical location of the sonic line in a meridional 
plane of the flow about a finite cone is of considerable interest 
for a theoretical study of the problem of axially symmetric 
flow. With this fact in mind, an interferometric analysis 
was made at several typical values of so as to determine 
the local Mach number fields about finite cones. Several 
examples of supersonic to subsonic shock-free compression 
are experimentally demonstrated. 

Experimental values of the shock-wave angle at the cone 
tip, particularly at values of where the flow field between 
the shock wave and the cone surface is transonic or subsonic 
in nature, are presented, and a comparison with the values 
from conical theory is shown. 

The conical solution for flow about a semi-infinite cone 
demonstrates that a conical solution does not exist if, for 
a given cone angle, decreases below a certain minimum 
This minimum is defined to be the for which 
shock-wave detachment occurs for a semi-infinite cone. 
Whether or not the shock-wave-detachment Mach number 
for a finite cone can be determined from conical theory is of 
considerable theoretical interest. Experimental values of 
the detachment distance of a shock wave from a finite cone 
tip, the distance obviously being zero at attachment, have 
been collected from several som'ces and the results analyzed 
in this report. 

The transonic equations of motion and boundary condi- 
tions as derived by Von Kdrmdn (ref. 9) for axially sym- 
metric flow require several assumptions as to the relative 
magnitude of various terms in the exact equations of motion 
and the related boundary conditions. To demonstrate that 
the transonic equations retain the important features of the 
exact equations, an approximate solution of the problem of 
conical flow about a semi-infinite cone has been developed 
employing the transonic equations and boundary conditions. 


A comparison of the exact Busemann solution and the ap- 
proximate transonic solution is presented in the report. 

The author of this report wishes to express his apprecia- 
tion for their helpful advice and criticism to Drs. H. W. 
Liepmann, J. D. Cole, and A. Roshko of the California 
Institute of Technology. The investigation was conducted 
under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 

SYMBOLS 

2k+\ 

a* velocity of sound for A/„ = 1.00 
b width of increments of region of integration 

Cj) cone drag coefficient; reference area is cone base area 
( z)* drag coefficient at Af« = 1.00 

Cp pressure coefficient, 

Cp* pressure coeflficient at M„ = 1.00 

c chord of cone 

d cone base diameter 

i, k integers 

I path length 

Ms cone-surface Mach number 
free-stream Mach number 

Mstc Mach number immediately downstream of a shock 
wave 

N number of outermost increment of region of integra- 
tion 

n index of refraction of air 

Po stagnation pressure 

p, surface static pressure 

p^ free-stream static pressure 

velocities in axial and radial directions, respectively 
r radial distance from axis of symmetry to point on 
light path 
cone base radius 

5 interferometric fringe shift 

u nondimensional transonic axial-velocity perturbation 

in appendix C; in appendix B 

u,=y.^ 

V nondimensional transonic radial-velocity perturba- 

tion in appendix C; in appendix A 
Vo V on surface of body 

w=r^ 

X axial distance downstream of cone tip 

y perpendicular distance from axis of symmetry to 

light path 

j8 shock-wave angle 

shock-wave angle at cone tip 
7 ratio of specific heats of air, 1.400 

6 axial distance from cone shoulder to shock wave 

7}=r/x 

r}s tangent of shock-wave angle 

B cone semiangle 

K Gladstone-Dale constant 

\o wavelength in vacuum of light employed 
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p density of air 

p(l) density of undetermined medium 

Pr density of reference medium 

p^ free-stream density 

(T=xjr 

<t> perturbation potential 

EXPERIMENTAL EQUIPMENT 

WIND TUNNEL 

The transonic wind tunnel at GALCIT is a continuous- 
flow wind tunnel. For supersonic testing, the test-section 
Mach number may be continuously varied over a wide range 
by altering the shape of one flexible w^ll. The test-section 
Mach number is varied by changing the area of a sonic 
throat downstream of the test section when subsonic tests 
are being performed. The test-section width is 4 inches and 
the height is 9 inches. The design of the flexible test-section 
wall is discussed in reference 16. 

INTERFEROMETER ' 

The interferometer used in the present investigation is of 
the Mach-Zehnder type. Both light paths of the inter- 
ferometer are passed through the wind-tunnel test section, 
one beam passing through the flow region under analysis 
and the other beam passing through the undisturbed flow 
in the test section upstream of the model. The fringe 
shifts due to density variations in the boundary layer are 
eliminated since both beams pass tlirough the boundary 
layer on the test-section walls and are affected almost 
equally. The growth of the boundary layer between the 
two beams is not compensated by the above arrangement, 
but the effect is of secondary importance. A detailed 
description of the GALCIT interferometer and a very com- 
plete bibliography on interferometer construction are given 
in reference 17. 

MODELS 

The models were conical-tipped brass cylinders of 20°, 
25°, and 30° semiangle. The base diameters were between 
0.30 and 0.50 inch. Thus, the Reynolds numbers for the 
tests, with the base diameter of the models as the reference 
dimension, varied from 55,000 to 143,000. The tips of the 
cones were made as sharp as possible and the maximum tip 
diameter of the dullest of the models was approximately }i 
percent of its base diameter. Also, the models were black- 
nickel plated to improve the photographic definition. The 
angle of attack and angle of yaw were adjusted to zero by 
equalizing the pressure on the cone surface at four annular 
points. 

EXPERIMENTAL TECHNIQUES 

INTERFEROMETRY 

An experimental investigation of flow phenomena is facil- 
ilitated by the employment of an interferometer to deter- 
mine the density fields in gaseous (or liquid) flows. The 
interferometer technique possesses the obvious advantage 
of eliminating the need for placing any type of probe into 
a flow region where the presence of the probe may com- 
pletely alter the undisturbed flow field. A disadvantage is 
also present, however, since the values of density are not 


immediately available as the test is in progress. A more 
serious disadvantage is the fact that the interferometer 
integrates the density values on its light paths (see appendix 
A) and, thus, the measurement of density is not localized 
but is influenced by inhomogeneities in the flow which may 
be well removed from the points of interest. 

Two general types of flow are amenable to interferometric 
analysis, namely, two-dimensional and axially symmetric 
flow. This paper is concerned solely with axially symmetric 
flow analysis. A discussion of the method employed to 
reduce the finite-fringe interferograms, such as figure 1, to 
density distributions may be found in appendix A. The 
method is essentially that outlined in reference 18. An 
excellent discussion of axially S 3 onmetric data reduction is 
given in reference 19 where several references to earlier 
papers in the field will be found. 







Figure 1. — Finite-fringe interferogram for 30“ semiangle cone cylinder. 

JW. = 1.280. 

These earlier papers are mainly concerned with evaluating 
the interferometer data-reduction techniques for axially 
symmetric flow by investigating the flow about cone cylin- 
ders at Mach numbers and cone angles where the Busemann 
conical solution was known to be valid. Reference 20 pre- 
sents some experimental results in the same general flow 
regimes as are investigated in this report. 

SONIC-LINE LOCATION BY WAVE REFLECTION 

The location of the sonic line in a meridional plane of an 
axially symmetric transonic flow can be determined experi- 
mentally by at least three distinct methods. The first 
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Figube 2. — Wave reflection from sonic line for 20° semiangle cone 
cylinder. A/ ^ = 1 . 297. ^ 

o Shock reflection doto, 



Figube 3. — Sonic-line location for 25° seiniangle cone cylinder. 

method is by static-pressure measurements, the second is 
by interferometric analysis, and the third is that of Mach 
wave reflection from the sonic line. It should be noted 


that at the point of reflection the Mach wave will be per- 
pendicular to the streamline direction through the sonic line. 

To locate the sonic line within the flow about a cone, a 
small-diameter probe was placed in the free-stream flow 
outside the cone shock wave. The probe was in a position 
such that the probe shock wave pierced the cone shock 
wave and entered the flow field about the cone. The shock 
waves formed by the probe closely approximate Mach waves 
at large distances from the probe and a typical wave reflec- 
tion is shown in figure 2. The perturbations in the flow 
about the cone caused by the waves do not appreciably 
affect the position of the sonic line as shown by figure 3 
where a comparison is made of the location of the sonic line 
as found by interferometric analysis and by the wave- 
reflection method. The probe method is much more con- 
venient than the pressure-measurement or interferometric 
method since the phenomenon may be observed with a 
schlieren system, so that the result is obtained visually. 

PRESSURE MEASUREMENTS 

Tlie pressure measurements in this investigation were 
made either on a micromanometer (accuracy of ±0.01 
millimeter of mercury) or on a nomograph Mach meter 
(ref. 21). 

EXPERIMENTAL RESULTS AND DISCUSSION 

GENERAL FLOW CHARACTERISTICS 

An analysis of the flow of a compressible fluid about an 
axially symmetric finite cone, that is, a cone cylinder, indi- 
cates that five distinct regimes of flow are possible. These 
regimes are given below. 

Regime I. — Regime I is subsonic flow at infinity with a 
region of locally supersonic flow downstream of the cone 
shoulder. A schlieren photograph of this type of flow is 
shown in figure 4. It should be noted that an extremely 
weak shock wave originates at the cone shoulder and 
terminates at the downstream “normal’’ shock wave. The 
forked appearance of the base of the terminating “normal” 
shock wave is an illusion caused by the axial symmetry of 
the flow. The light rays which pass near the surface of 
the body in the region of the rearward branch of the “fork” 
also pass through the outer portion of the shock and the 
spurious rearward branch is caused by the light-ray deflec- 
tions in the outer portion of the shock wave. A meridional 
section of the shock wave actually includes only the front 
branch of the fork. 

Regime n. — Regime II is supersonic flow at infinity 
with a detached shock wave and subsonic flow between the 
shock wave and the cone. Figure 1 is a finite-fringe inter- 
ferogram of this type of flow. 

Regime HI. — Regime III is supersonic flow at infinity 
with an attached curved shock wave and subsonic flow 
between the initial portion of the shock wave and the cone. 
A schlieren photograph of this flow is shown in figure 5 (a). 
Taylor and Maccoll’s original paper on conical flow (ref. 13) 
includes a schlieren photograph of an attached curved shock. 

Regime IV. — Regime IV is supersonic flow at infinity with 
an attached shock wave and mixed supersonic and subsonic 
flow between the shock wave and the cone. A schlieren 
photograph of this flow is shown in figure 5 (b). 
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Figure 4. — 20° semiangle cone cylinder. = 0.942. 


Regime V. — Regime V is supersonic flow at infinity with 
completely supersonic flow between the attached shock 
wave and the cone surface. The Busemann solution applies 
in this regime and has been verified experimentally in refer- 
ences 13, 14, 19, and 22. 

LOCAL MACH NUMBER CONTOURS 

The local Mach number contours in a meridional plane 
for the flow about a 25° semiangle cone are shown in figure 
6 for flow regimes II, III, and IV. The local Mach number 
contours for a 30° semiangle cone in regime II are shown in 
figure 7. These data were obtained by interferometric 
analysis as discussed in appendix A. 

The local Mach number contours should be normal to the 
cone surface since the cone surface is a flat boundary and 
any pressure gradient at the surface must be parallel to the 
flat surface. However, near the shoulder of the cone 
cylinder the surface is curved by the effect of the corner 
expansion on the boundary layer, and thus the local Mach 
number contours are not quite perpendicular to the cone 
surface immediately ahead of the shoulder. 

SONIC-LINE LOCATION 

The location of the sonic line can be determined by inter- 
ferometric analysis, but a more useful method, in the present 
investigation, was the wave-reflection method. The loca- 
tion of the sonic line in regimes II, III, and IV is discussed 
below. 


Regime II. — Examples of the sonic-line location in regime 
II are shown in figures 6 (a) and 6 (b). The sonic line 
originates at the cone shoulder and terminates on the 
detached shock wave. In figure 7 it may be seen that a 
region of supersonic to subsonic compression exists on the 
outer portion of the sonic line. The sonic line actually 
originates slightly upstream of the cope shoulder. This 
effect is due to the rounding of the cone shoulder by the 
surface boundary layer. 

Regime III. — Figure 6 (c) illustrates the case of the flow 
with a nearly attached curved shock wave. Again, a small 
region of supersonic to subsonic compression is present on 
the outer portion of the sonic line. The free-stream Mach 
number is slightly less than the detachment predicted 
by the exact conical theory. The question of experimental 
detachment Mach number is discussed subsequently. 

Regime IV. — Several examples of the sonic-line location 
in regime IV are shown in figures 3, 6 (d), and 8. Figure 3 
shows the location as determined by interferometric means 
and as determined by wave reflection. The sonic line again 
originates at the corner and now terminates at the cone tip 
and not on the shock wave as in regimes II and III. A 
shock-free supersonic to subsonic compression occurs on the 
forward portion of the sonic line. The location of the sonic 
line for a 20° semiangle cone is shown in figure 8. The 
agreement between the theoretical and experimental location 



(a) A/eo = 1.328. 

Figure 5. — 25° semiangle cone cylinder. 
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is satisfactory near the tip of the cone and for some distance 
downstream of the tip.^ 

The question of smooth shock-free supersonic to subsonic 
compression has been the subject of much discussion in 
recent years. The above experimental results demonstrate 
that such a flow is possible. However, the smooth com- 
pression is not of primary importance, but rather the 
conditions under which it occurs. These conditions are that 
the sonic surface bounds a zone of subsonic flow completely 
enclosed by a region of supersonic flow and a solid surface. 


/ 



(b) M„ = 1.415. 
Figure 5. —Concluded. 


As an example of non-shock-free supersonic to subsonic, 
compression, consider the flow past a two-dimensional air- 
foil at high subsonic speeds. The local supersonic zone on 
the airfoil is terminated by a shock wave and smooth com- 
pression through sonic velocity does not occur. In the two- 
dimensional case, however, the supersonic zone is bounded 
by a subsonic region and a solid surface. This is the op- 
posite arrangement to that in the flow about a cone, in 
regime IV, where shock-free supersonic to subsonic com- 
pression does occur. 

The above considerations illustrate that the existence (or 
stability) of shock-free compression through sonic velocity 
may not be a local plienoiiienon but may depend on the 
arrangement of the complete flow field. 

* The existence of this type of flow is indicated by the experimental results of Taylor 
and Maccoll (refs. 13 and 14) and \vas also discussed by Tsien (ref. 23). 



(a) A/„ = 1.229. 

Figure 6. — Mach luuuber field for 25° semiaugle cone cylinder. 



Figure 6 . — Continued. 
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(c) = 1.273. 

Figure 6.— Continued. 



Figure 7.— Mach number field for 30® semiangle cone cylinder. 





(d) = 1.401. 

Figure 6 . — Concluded. 



. Figure 8.— Sonic-lino location by wave renection for 20° semiangic 

cone cylinder. 


.{50080—50 2 
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SHOCK-WAVE ANGLE 

The angle of the attached shock wave at the nose of the | 
cone was determined for a 20° and a 25° semiangle cone in | 
flow regimes III, IV, and V. The values are shown in figure ! 
9. Similar experiments are reported in references 13 and 14. 
Reference 22 presents data for one cone angle at one Mach 
number in regime III and one Mach number in regime IV. | 
The agreement between the exact theory and the experi- 
mental values at the, cone tip is very good even in regimes 
III and IV where the exact theory is not applicable for the 
complete finite cone. 



SURFACE MACH NUMBER DISTRIBUTION 

The distribution of the surface Mach number M, on a 25° 
semiangle cone for various values of M* is shown in figure 
10 and that on a 20° semiangle cone, in figure 11. Several 
characteristics of these distributions are of particular interest. 

(1) Surface Mach number near the shoulder deviates from 
the Busemann conical values as soon as M^=1.00 is attained. 
Surface Mach number near the cone tip agrees quite well 
with the conical values until the theoretical detachment Moo 
occurs. At the corner M, should, except for boundary-layer 
effects, alwa\^s be sonic if M» forward of the shoulder is 
subsonic. 

(2) As Moo approaches 1.00 from the subsonic or from^tlie 
supersonic regimes, M, at a particular chordwise station ap- 
proaches a constant value. This behavior implies that 

=0 



Figure 10. —Surface Mach number on 25® semiangle cone cylinder. 



Figure 11.— Surface Mach number on 20° semiangle cone cylinder. 


The same behavior of M, on two-dimensional sections was 
noted in reference 24, and thus the concept of stationary 
values of Ms at M« = 1 .00 is established for two-dimensional 
and axially symmetric flow. Since these two cases represent 
limiting cases of the flow about general three-dimensional 
bodies, the stationary M, concept can probably be applied 
quite generally if suitable care is taken in choosing the range 
of Moo about Moo = 1.00 in which the so-called '‘M, freeze” 
is applicable. 

(3) As Moo progresses from a subsonic value through 
Moo = 1.00 and on to a value in regime V, the M, at a par- 
ticular chordwise station probably varies quite smoothly 
with no abrupt variations, even at attachment of the shock 
wave, except for a region quite near the tip where large 
variations may occur when the shock wave attaches. 


DRAG COEFFICIENTS 


The drag coefficients for the 20° and 25° semiangle cones 
are shown in figure 12. The values at Moo = 1 00 were de- 
termined by extrapolating the M, data in figures 10 and 11 
to Moo = 1.00. 

Using the concept of stationary values of M, at Moo = 1.00, 
the drag-curve slope at Moo = 1.00 becomes (see ref. 4 and 
appendix B) 


dCu 

dM^ 


_ ^ * 

T+1 7+1 ^ 


dM^ 


A/oo^l.OO 


A/oo=1.00 


( 1 ) 
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where is tlie drag coefficient at Ma, = 1.00. The first 
term 4/(7+!) of the drag-curve slope is derived from the 
first-order term of the pressure coefficient while the term 
[2/(7+ !)](?/)* represents the contribution of the second-order 
terms. The magnitude of the second-order term [2/(7+ 
1)]^D* is shown by the difFerence in slope of the pairs of 
lines drawn through Co* in figure 12. 

The experimental results also indicate that 


clMJ 


=0 

A/ CO ~ 1. 00 


This then implies (see appendix B) that 


cPCo 


dU^ 


A/00 =1.00 


127+4 107+6 ^ * 

(7+1)^ 7t+1)^ 


and an estimation can then be made of the range about 
Moo = 1.00 where equation (1) is valid. 


SHOCK-WAVE DETACHMENT 

Conical-flow theory indicates that for a given cone angle 
of a semi-infinite cone a certain minimum Mr^ is reached 
below which a conical solution is no longer possible. This 
value of Moo is defined to be the shock-wave detachment Moo. 
However, a finite cone introduces a characteristic length into 
the problem so that curved attached shock waves, which 
would provide the necessary pressure gradient to turn the 
flow near the cone tip, may exist at values of Mco less than 
the conical detachment Mc^.* 

Present experimental results indicate only that shock-wave 
detachment for a given cone angle does not occur at an 
Moo greater than that predicted by conical theory. A col- 
lection of data from references 22, 25, 26, and 15 is shown in 
figure 13. The ratio 6/rf,^where 5 is the center-line distance 
from the shock wave to the plane of the cone shoulder and 
d is the body diameter at the shoulder, that is, the sonic 
point, is seen to approach asymptotically the value of hfd 
at attachment. The asymptotic behavior of hjd complicates 
the fairing of the proper curve of hjd versus particularly 
in view of the paucity of experimental points in the immediate 
vicinity of shock-wave attachment. 

In reference 27 data are presented for the shock-wave 
detachment distance of several cone angles at Mo, =2.45. 
A discrepancy was found between the experimental and 
theoretical values of the cone angle at which shock-wave 
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Figure 13. -Shock-wave detachment distance for finite cones. Data 
are primarily from references 15, 22, 25, and 26. 

detachment occurs for a fixed value of M^, detachment ap- 
pearing to occur at a cone angle slightly greater than that 
predicted by the conical theory. This behavior would cor- 
respond to shock- wave detachment for a fixed cone angle 
occurring at a value of M« less, than the theoretical conical 
value of M«. Again, however, the discrepancy may be 
caused by the manner in which the experimental curve was 
faired. 

Thus, the experimental results appear to indicate only that 
shock-wave detachment for a finite cone occurs at a value of 
Mco either less than or equal to the value of M^ predicted by 
conical theory but not at a larger Mco. The fact that shock- 
wave detachment does not appear to occur at a value of 
greater than that predicted by the conical theory indicates 
that the presence of a boundary layer on the cone tip does 
not affect the conicity of the flow near the cone tip to the 
extent of precipitating detachment of the shock wave. 

Figure 13 also indicates that when the shock-wave detach- 
ment distance is large, the position of the shock wave is 
dependent only on the diameter of the cone at the shoulder 
or sonic point and is independent of the cone angle. When 
the shock wave is quite near the cone tip, however, the de- 
tachment distance is also dependent on the cone a,ngle. This 
manner of behavior of the shock-wave separation distance 
was discussed by Busemann (ref. 28) and was shown experi- 
mentally for two-dimensional wedge sections by Griffith 
(ref. 29). 

TRANSONIC SIMILARITY 

The transonic-similarity rules for the drag coefficient and 
pressure coefficient, as derived in reference 10, cannot be 
checked by the experimental results of this report. The 
derivation assumes that the cone-surface boundary condition 
is the approximate tangency condition which is valid for 
relatively small angles. A 20"^ semiangle cone is the mini- 
mum-angle cone for which detached shock-wave flow can be 
obtained in the transonic wind tunnel, and thus the experi- 
mental models were 20°, 25°, and 30° semiangle cones. The 
experimental cone angles are much larger than the cone angles 
for which the approxunate tangency condition is reasonable, 
and, therefore, the transonic-similarity rules of reference 10 
are not applicable. 
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THEORETICAL CONSIDERATIONS 

At the present time, theoretical solutions have not been 
found to describe the flow about a finite cone for the com- 
plete Mach number range. Theoretical solutions are 
available for only two Mach number regimes, namely : 

(1) Exact conical theory may be applied if the surface 
Mach number is greater than srmic 

(2) At Moo = 1.00, Yoshihara (ref. 11) lias calculated 
the flosv about small-angle cones by a relaxation 
technique 

No solution has been determined if Moo is less than sonic. 
However, Van Dyke's second-order supersonic theory and 
technique (ref. 30) possibly can be applied to the subsonic 
case since, if the appropriate changes of sign are made in 
the particular solution found by Van Dyke for the super- 
sonic case, the particular solution becomes valid for the 
subsonic case. 

A solution remains to be found for the regime between 
M« = 1.00 and the value of M* where Ms becomes equal to 
LOO. The problem would be greatly simplified if the 
transonic equations could be employed. To test the feasi- 
bility of the approximations inherent in the transonic equa- 
tions, an approximate solution has been found for conical 
flow about a semi-infinite cone using the transonic equations. 
The details of the solution are presented in appendix C. 

The solution is compared with the exact conical theory in 
figures 14 and 15. Figure 14 shows the comparison between 
the shock-wave angles predicted by the exact theory and by 
the transonic approximation. The surface Mach number 
comparison is shown in figure 15. From figure 15 it can be 
seen that the transonic approximation is qjLiite satisfactory 
and is probably better than slender-body cone theory, since 
slender-body cone theory does not consider the presence of 
the conical shock wave. Also, figures 14 and 15 show the 
agreement of the cone angle at shock-wave detachment as 
found from the transonic solution and from the exact theory. 

The above comparison of the exact conical solution and 
the approximate solution indicates that the transonic equa- 
tions contain all the terms of importance in the exact 
equations for conical flow about cones, so that the transonic 
equations may be employed with confidence in the range of 
Moa from M„ = 1.00 to an M* for which Ms^l.OO. 

SUMMARY OF RESULTS 

The following results were obtained from an investigation 
of transonic flow past cone cylinders: 

1. The variation of drag coefficient Co with free-stream 
Mach number M* was determined experimentally. The 
slope of Cd versus Mo, at M«, = 1.00 agrees with the theo- 
retical prediction. The deviation of Cd versus M^ from the 
conical flow value of Cd when M,<1.00 is demonstrated. 

2. The experimental results for the shock-wave angle, 
sonic-line location, and surface Mach number in the region 
near the cone tip indicate that the flow is conical near the 
tip of a finite cone even when the surface Mach number is 
less than sonic. The surface Mach numbers for the rest of 
the cone deviate from the exact conical values when Ms< 1 .00. 
Also, a case of shock-free supersonic to subsonic compression 
is demonstrated experimentally. 



Figure 14. — Shock-wave angle for semi-infinite cones. 
u„=0.25; A/„ = 1.12. 
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Figure 15. — Transonic conical solution; surface Mach number versus 
cone semiangle. w„ = 0.25; A/„ = 1.12. (A/„ is defined as M^ = 


3. An approximate solution for transonic conical flow has 
been developed and the agreement with the exact conical 
theory indicates that the axially symmetric transonic equa- 
tions retain the important features of the exact equations. 

4. Present experimental values of the detachment distance 
of a shock wave from a finite cone tip do not demonstrate 
agreement with the detachment Mach number predicted by 
conical theory for a semi-infinite cone and the question of 
shock-wave detachment from a finite cone remains undecided. 
California Institute of Technology, 

Pasadena, Calif., A'pril 16, 196S. 


APPENDIX A 


REDUCTION OF INTERFEROMETER DATA 


An interferometer determines the advancement or retarda- 
tion of a light wave in a medium with respect to a coherent 
light wave in a reference medium. Since the wave velocities 
are a function of the indices of refraction of the respective 
mediums and consequently of the densities of the mediums, 
it may be shown that 

(Al) 

where 

p(/) density of undetermined medium 

Pff density of reference medium 

I path lengtli 

K Gladstone-Dale constant 

Xo wavelength in vacuum of light employed 

S fringe shift; in case of finite-fringe interferograms 
this is ratio of displacement of a fringe to interval 
between undisturbed fringes 

In equation (Al) it has been assumed that the light beams 
traverse identical geometrical paths, so that refraction, if 
present, is neglected. Also, the relationship between the 
index of refraction n and density is assumed to be 

71=1 -{-Kp (A2) 


If 7i=l + a where a<^l, equation (A2) is obtained by 
linearizing the Lorentz law. 

For the axially symmetric case, the fringe shift for a light 
path perpendicular to the axis of symmetry becomes 



[p{r) — pR]rdr 



where r is the radial distance from the axis of symmetry to 
a point on the light path and y is the perpendicular distance 
from the axis of symmetry to the light path. 

In the present investigation, the density field was bounded 
by a shock wave at a distance 7/, from the axis and the refer- 
ence density was the free-stream density p«, ; thus. 



[p(r)— poo]rrfr 

(^2_y2)l/2 


(A3j 


Weyl (ref. 31) introduced the assumption that S(y) is a 
linear function of in a small interval of y. The validity 
of this assumption for the present investigation is indicated 
by the parabolic nature of the typical fringe-shift curves 
shown in figure 16. If the substitutions 


71 = / 


Figure 16. — Interferogram fringe shifts for 25° semiangle 
cone cylinder. 


are made, equation (A3) becomes 

S{u)=— [piv) — p<x>]dv 

— u 


(A4) 


This is analogous to the solution of Abeks problem. The 
solution of equation (A4) for p is 


jdSl^n 

^ ^lu — w 


(A5) 


where w=r^. A complete proof of the solution may be 
found in reference 19. 

If the region of integration in equation (A5) is divided into 
equal increments in r of width 6, 


Us=yJ^ 


Ti=ib 
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where i is an integer. A numerical evaluation of equation 
(A5) is then 



pW_1+ 

(A6) 


Poo TTKCLpao k^i 


wliere 

r,=Nb 


and 




2jt+l 


The above solution of the problem is essentially that of 
reference 18 and a table of A^i for 50 intervals will be found 
in reference 18. 

From the density ratio determined by equation (A6) the 
local Mach number may be computed. An approximate 
correction to the local stagnation pressure po' downstream of 
the shock wave was made by assuming that on the cone 
surface po' was determined by the nose shock-wave angle 
and at a given chordwise station po' varied linearly with the 
value of the shock wave at the given station. If the approx- 
imate streamlines are calculated, as in reference 20, a more 
refined correction is obtained. From the experimental 
values of the local Mach number a topographic map was 
plotted, and from this map the desired Mach number 
contours are found. 

CONICAL FIELDS 

If it is desired to determine whether or not a given axially 
symmetric field is conical, a simple test can be made. Re- 
turning to equation (A3) it is assumed that 

p=p(v) 

where rj^rlx and x is the axial distance from the conical 
origin. Then equation (A3) becomes 

S{y)_2K p* [p{r))—Pa ,] ydr) 


where Vs is the tangent of the shock-wave angle and vi=ylx- 
Thus, if the field is truly conical 

and a plot of S(y)/x versus yjx for various values of x will 
yield a group of coincident curves. Examples of this tech- 
nique are shown in references 19 and 20. It is interesting to 
note that in reference 20 figure 6 (b) indicates conical flow 
near a cone tip for flow regime IV, that is, a 35° scmiangle 
cone at Moo = 1.87, using the above technique. 

SUBSIDIARY CONSIDERATIONS 

Model size. — From equation (x\3) it is evident that the 
fringe shift at a particular chordwise and radial station is a 
linear function of the model size for fixed values of density. 
Thus the model should be as large as is compatible with the 
test-section dimensions with regard to blocking and so forth. 

Finite fringe spacing. — The fringe spacing in the undis- 
turbed field must be such that a sufficient number of data 
points may be determined between the shock wave and the 
cone surface at a particular chordwise station. However, 
for a given fringe shift S the displacement of the fringe is 
proportional to the undisturbed fringe spacing, and the 
accuracy of the fringe data will be improved by increasing 
the undisturbed spacing. A compromise must be effected 
between the desire for many fringe shift points at a given 
chordwise station and the accuracy of the individual points. 
In the present investigation, this compromise precluded a 
study by interferometry of the flow properties in the immedi- 
ate vicinity of the cone tip when the shock wave was attached. 

Accuracy. — The accuracy of the interferometric method is 
affected by refraction, inhomogeneities in the reference flow, 
the numerical approximation, and so forth. An estimate of 
the accuracy can be obtained by noting the comparison of 
interferometer data and shock-reflection data in figure 3 
and the values of local Mach number behind the shock wave 
indicated in figures 6 and 7. 

Comprehensive discussions of the accuracy of the method 
may be found in references 18 and 19. 


APPENDIX B 

VARIATION OF PRESSURE AND DRAG COEFFICIENTS NEAR Mco=1.00 


The pressure coefficient 



where Ps is the surface static pressure and is the free- 
stream static pressure may, in transonic flow, be written as 



where M, is the surface Mach number. In equation (Bl) it 


has been assumed that 

Vo 

that is, the stagnation pressure loss across any shock waves 
may be neglected. 

Defining f and $ as 

?=M,2 

equation (Bl) becomes 
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For a fixed body geometry, M, is a function of only; 
therefore, 




Thus, the first a?id second total derivatives of Cp with respect 
to f are 


c/f 


m 


Tlic partial derivatives /ff and f( may be calculated from 
equation (B 2 ) and thus 


\dM^J 7+i 



127+4 107+6 ^ * 

dMjJ (7+1)^'' (7+ If 


(B5) 

(B 6 ) 


and 


df 


h+Vn |+./« (^) +/« 




(B4) 


The fore drag coefficient of a finite cone is given by 



In general, and d^^/df^ are not known; however, an 

argument presented by Liepmann and Bryson (ref. 24) 
shows that, for ^'=1.00, d^/df=0. The same argument, 
namely, that has a stationary value at Mco = 1 . 00 , 
cannot be used to evaluate However, an inspection 

of the experimental data (see figs. 10 and 11) indicates that 
(he curve of vei-sus video lias an inflection point at Meo = 
1.00 and, thus, d^f/df^=0 at f=1.00. 

Using equatio/is (B3) and (B4) and the above argument, 
the derivatives of Cp become 



where ( )* indicates evaluation at f= 1 . 00 , that is, at M«,= 

1 . 00 . The derivatives of Cp with respect to are then 


where oL={rlrtY^ r is the local cone radius, and ri, is the cone 
base radius. Using equations (B5) and (B 6 ), the derivatives 
of Cd are 


and 


/ df^ J 2_ 

\dM^) 7+1 7+1 


(B7j 


/d^C^Y__127+4 107yf6 ,, , 

\dMj) (7+l)2^(7+l)2^^^ 


(B 8 ) 


A comparison of equation (B7) and experimental data is 
shown in figure 12 . Equations (B 5 ) and (B 7 ) were given 
previously by Bryson (ref. 4). 

It should be noted that the first term 4 /( 7 +!) of equation 
(B5) may be derived from the linearized transonic Cp (ref. 1 ) 
which is 


Cp^2 


MJ-MY 

(7+1) 


(B9) 



The second term — [ 2 /( 7 + 1 )]( 7 ^* of equation (B 5 ) is then of 
the nature of a second-order correction term. However, if 
{(PCpIdMcY)"^ is computed from equation (B 9 ), the result 
does not agree with the first term of equation (B 6 ) and is, 
in fact, of opposite sign. 


APPENDIX C 


TRANSONIC APPROXIMATION FOR CONICAL FLOW 


In discussing axially symmetric transonic flow, the follow- 
ing approximations to the exact equations and boundary 
conditions are employed. 

If Qx is the velocity in the axial or a:-direction and qr is the 
velocity in the radial or r-direction it is assumed that 

qx=^a*+<t>x 


the continuity equation is approximated by 

UUx—Vr—-V=^ (Cl) 

r 

and the condition of irrotationality becomes 

Ur-^Vz=0 (C 2 ) 


<lT=^r 


where a* is the velocity of sound at Mco=1.00, 0 is the 
perturbation potential, and (t)r<^a*. Then, defining 


u- 


^ 7+1 

a* 

_ 7+l 

a* 


<t>x 


<Ar 


The exact tangency condition on the body surface is 
replaced by 

«?o=(7+l) tan 0 (C3) 

where Vo is v on the surface of the body and 6 is the inclination 
angle of the body surface. The above relations are derived 
in greater detail in reference 9 . 

To test the feasibility of the approximations inherent in 
equations (Cl), (C 2 ), and (C3) an approximate solution for 
the flow about a semi-infinite cone will be developed. This 
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approximation will then be compared with the Busemann 
solution of the exact equations. 

CONICAL SOLUTION 

Assuming that u and v are functions of <t where 

X 

<T== — 

r 


equations (Cl) and (C2) become 

uu«r+oi;<r— r=0 

A solution of the form 

u=J{v) 

will be sought. From equation (C4) 


(Cla) 

(C2a) 

(C4) 


dcr dv da 


but from equation (C2a) 


dv 


du 


""da 


Therefore, 


du 

dv 


(C5) 


The relationship in equation (C5) is exact (see ref. 32). 
Differentiating equation (C5) with respect to a yields 


dv 


1 


da 2 d^u 


and, therefore, 


du 


1 


da 


dv^ 


(C6) 


(C7) 


Substituting equations (C5), (C6), and (C7) into equation 
(Cla), the result is 


(Pu , du /duV 


(C8) 


The following conditions are imposed on equation (C8). 

(a) At the shock wave : If Us and Vs are the values of u and v 
at the shock wave, then, from the transonic-shock polar 
relation, 

(C9) 




where 


u„=MJ-l 
The shock-wave angle is given by 

Vs 


cot )3= 


u—u, 


hence, from equation (C5), 

du 
dv 


u„—Us 


(CIO) 


(b) On the cone surface : The tangency condition must be 
satisfied; hence, 

«^o=("f+l) tan d (Cll) 

where Vo is v on the cone surface and B is the cone semiangle. 
From equation (C5) 

du 


dv 


= — tan B 


(C12) 


An exact solution of equation (C8) has not been found. 
However, if the right-hand side of equation (C8) is assumed 
to be small, an iteration solution may be found. 


FIRST APPROXIMATION 


As a first approximation to the solution, set the right-hand 
side of equation (C8) equal to zero; that is. 


d ( du\^ 
dv \ dv/ 


The solution of equation (Cl 3) is 

u^A loge Bv 


(C13) 


(C14) 


Applying conditions (C9) and (CIO) to this solution, equation 
(Cl 4) becomes 


u,) loge “ 


SECOND APPROXIMATION 


(C15) 


As a second approximation to the solution of equation 
(C8), the right-hand side of equation (C8) is assumed to be 
given with sufficient accuracy by equation (Cl5). Hence, 


d ( du\ 
dv \ dv / 


Usiu^—UsY , {u^—UsY 1 ^^ V 


(C16) 


The solution of equation (Cl 6) is 

M = — M.) + (W,—U.) loge ^J+ClogeZ?t> 

(C17) 

Applying conditions (C9) and (CIO) to this solution, equation 
(Cl 7) becomes 


u^ 




5«,) 


(u^+5u,)(u„—u,) , V u^(u 


(C18) 


The values of u, and Vs appearing in equation (Cl8) are 
not independent of the cone semiangle B. Applying condition 
(Cl2) to equation (Cl8) and solving for Vsj the result is 


loge Vs=loge (y+1) tan B 

( 

\u^—u, 2/ 


[(7+1) tan eY (it^+5u,) 
(u„+u.) (u„— u,)® 


2(t+1)^ ta.n*e 
{u„-u,y 


(C19) 
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The transonic-shock polar relation also must be satisfied; 
hence, 

= (C20) 

Because of the nature of equation (Cl9), explicit solutions of 
equations (Cl 9) and (C20) for and Vs in terms of and S 
have not been found. However, a solution may be found 
graphically. 

The values of Us and Vg tlius determined for a given value 
of 6 and may be introduced into equation (Cl 8). If the 
value of u on tlie surface Uo is desired, then substitution of 

i^=ro=(7+l) tan d 


The usual transonic approximation reduces equation (C21) to 

Mg^==l+Uo (C22) 

Equation (C21) is shown in figure 15. The agreement with 
the exact theory is quite satisfactory. 

Apple curve. — Two sets of values of Ug and Vg will satisfy 
equations (Cl 9) and (C20). The two solutions correspond 
to the “strong’’ and “weak” shock-wave solutions predicted 
by Busemann. If the values of Uo and Vo for all possible 
cone angles and for both types of solution are plotted in the 
hodograpb, the resultant curve is Busemann’s “apple” curve. 
The apple curve found in the second approximation is shown 
in figure 17. 


in equation (CT8) yields Uo. 

It should be noted that, for a given value of d and Ua,y two 
sets of values for Ug and Vg are found. These correspond to 
the “strong” and “weak” shock waves, 

COMPARISON OF SECOND APPROXIMATION AND BUSEMANN CONICAL 

SOLUTION 

Wave angle. — The wave angle is determined by the 
values of Ug and Vg since 


The degree of agreement is apparent in figure 14. 

Surface Mach number. — In terms of Uo and Vo, the surface 
Mach number is 



where 


(C21) 





second approximation. 
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